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Abstract 



Unconventional way of handling the perturbative series is presented with the 
help of Heaviside transformation with respect to the mass. We apply Heaviside 
transform to the effective potential in the massive Gross-Neveu model and carry out 
perturbative approximation of the massless potential by dealing with the resulting 
Heaviside function. We find that accurate values of the dynamical mass can be 
obtained from the Heaviside function already at finite orders where just the several 
of diagrams are incorporated. We prove that our approximants converges to the 
exact massless potential in the infinite order. Small mass expansion of the effective 
potential can be also obtained in our approach. 
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1 Introduction 

Even if the proof of dynamical massless symmetry breaking requires genuine non- 
perturbative approaches, it does not necessarily mean that the perturbative expan- 
sion is totally useless. There is the possibility that non-perturbative quantities in 
the massless limit may be approximately calculated via perturbative approach. The 
purpose of this paper is to explore the possibility and show a concrete affirmative 
result by re- visiting the Gross- Neveu modeP. 

Let us consider the effective potential of the Gross- Neveu model. As is well known, 
ordinary massless perturbation expansion gives infra-red divergences and to cure the 
problem one must sum up all the one-loop diagrams. Then the summed result reveals 
the non-trivial vacuum configuration of < > and the dynamical generation of 
the mass. 

The point we like to note is whether such a non-perturbative effect needs, in the 
approximate evaluation, the infinite sum of perturbative contributions. To resolve 
the issue, we deal with a truncated series Vpert, without conventional loop summation, 
and study the approximate calculation of the effective potential at m = 0. 

A naive way of approximation would go as the following: To get around the 
infrared singularity we turn to the massive case and probe Vpertio'j'm) at small m. 
Since the limit, m ^ 0, cannot be taken in Vperticjin), we may choose some non- 
zero m {— m*) and approximate the effective potential V{a,'m — 0) by Vpert{c,m*). 
However, the problem is that Vpert{o', m) is not valid for small enough m. This is the 
place where the Heaviside function comes in. Our suggestion to resolve the problem 
is to contact the Heaviside transformation of V{a, m) with respect to the mass^'^. 

Heaviside transform of the effective potential, V , is a function of a and x which is 
conjugate with m. Then, the key relation is that lim^_»o V{a, m) — Imix^^ ^^(c, x). 
Of course this is valid only when the both limits exist and do not apply for Vpe.rt a-nd 
its Heaviside function, Vpert-, because those functions diverge in the limits. However 
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there arises the possibihty that V{o-, oo) and hence V{a, 0) may be well approximated 
by putting some finite value of x into Vpert- This is because Vpert has the convergence 
radius much larger than that of Vpert- Although Vpert shares the similar infra-red 
problems with Vpert, we will find that Vpert is much more convenient in this kind 
of massless approximation. Actually we will demonstrate that, at finite perturba- 
tive orders where just the several of Feynman diagrams are taken into account, the 
accurate dynamical mass is obtained via the Heaviside transform approach. 

Throughout this paper, we use dimensional regularization''. We confine ourselves 
with the leading order of large N expansion and N is omitted for the sake of sim- 
plicity. 

2. Heaviside transform with respect to the mass 

In this section we summarize basic features of the Heaviside transform and illus- 
trate our strategy by taking a simple example. 

Let 0(m) be a given function of the mass m. The Heaviside transform of 0(m) 
is given by the Bromwich integral, 

^(^) = / ^ ^(H, (1) 

where the vertical straight contour should lie in the right of all the possible poles 
and the cut of Q,{m)/m (In (1), the real parameter s specifies the location of the 
contour). Since Q{m)/m is analytic in the domain, Re{m) > s, Q{x) is zero when 
x < 0. It is known that the Laplace transformation (of the second kind) gives the 
original function as, 

/oo 
dxexp(—mx)d{x). (2) 
-oo 

Since Q(x) — ior x < 0, the region of the integration effectively reduces to [0, oo). 
It is easy to derive the relation, 

lim fl(m) = lim fl(x), (3) 
3 



where the both hmits are assumed to exist. As noted before, the point of our scheme 
consists in utihzing Cl to approximate the massless value of il, fl{0), by relying upon 
(3). 

To illustrate our strategy based on (3), let us consider a simple example. Given a 
following truncated series in 1/m, 

him) = E (4) 

n=0 

we try to approximate the value of /(m) = fooijn) = (1 + rn)~^ at m = 0, /(O) = 1, 
by using information just contained in the truncated series (4). Since the convergence 
radius, p, of fooiin) is unity, we cannot have approximation better than 1/2 from 
fLijn)- However, the state changes if we deal with its Heaviside function. 
The Heaviside transform of fLijn) is given by 

where 

Prom (5) it is easy to find that f{x) — {l — e~^)9{x) and (3) holds for / and /. For our 
purpose it is crucial that p = oo for f^o while p = 1 for /oo- The infinite convergence 
radius ensures us to probe the large x behavior of / by Jl to arbitrary precision by 
increasing perturbative order. Due to the truncation, however, diverges as x — > oo. 
Then, in approximating f{oo) and therefore /(O), we stop taking the limit and input 
some finite value into x. The input value of x, say x*, should be taken as large as 
possible in the reliable perturbative region in x. At this place we understand that 
the good convergence property of is one of the advantage of Heaviside function. 

Since the upper limit of perturbative region is not a rigorously defined concept, 
we determine the input value x* in heuristic way. Our suggestion to fix x* is as the 
following: The series (5) is valid for small x but breaks down at large x. The breaking 
appears as the domination of the highest term in Jl which leads to the unlimited 
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growth or decreasing of the function (see Fig.l). Thus x* is located somewhere around 
the beginning of the dominating behavior. Then, for odd L and large even L, we 
find the plateau region just before the domination and that the region represents the 
end of the perturbative regime. Thus, we choose the stationary point in the plateau 
region as representing the typical violation of the perturbation expansion. Hence we 
fix X* by the stationarity condition, 

dhix*) 



dx* 



0. (7) 



The condition (7) reads as 



and reduces to 

EH^ = 0. (9) 

n=0 

The solution exists for odd L and it varies with L. We find from (9) that the solution 
X* tends to oo as L — > oo. More precisely the solution scales for large L as 

X* ^ h. (10) 

We have explicitly done the numerical experiment to several higher orders and 
obtained the result for L = 1, 5, 9, 13, 17, 

hix*) = 0.5, 0.850675, 0.953301, 0.985166, 0.995251, (11) 

at X* = 1,2.18061,3.33355,4.47541,5.6112, respectively Thus, the sequence gives 
good approximation of the exact value. In this toy model, one finds that fhi^*) 
converges in the L — > oo limit by using the scaling relation (10). 

Up to now we have concentrated on approximating the massless value. We here 
point out that our scheme is capable of constructing the small m expansion of Q(m), 
that is, the scheme allows the approximation of the function itself when m is small. 
Consider in general the approximation of the derivatives at m = 0, 



(12) 

m=0 



which is needed when one constructs the small m expansion of fl{m), 

Q(m) = Q(0) + ^^P-m + ^^-^m' + ■■■. (13) 

The coefficients, O('=)(0) {k = 1, 2, 3, • • •), can be approximated as follows. The start- 
ing formula is that 

0('=)(m)^ £jti-tr^i^l (14) 



where we used 



dm ox m ' 



Here Ti. above the arrow represents the Heaviside transformation. Hence, from the 
agreement condition (3) we find 

lim nW(r.) = Ita = £ (16) 

Now in our perturbative approach we use (the truncated series at the order L) 
for real Q. Then, we show that we can simulate ak{oo) by Q;fc(x^) where x], may be 
fixed following the same reasoning we presented for the case of /(O) approximation. 
Namely we guess that the break down of the perturbative expansion is represented 
by the plateau region, if it exists, just before the unlimited growth of the size of the 
function. Therefore we use stationarity condition which reads. 



doik 



dx 



= 0' (17) 



and find that xl satisfies the same condition as that for x*. Thus the solution of 
(17) is universal for all k and fixes the coefficients of the small m expansion to all 
orders. This is desirable since the uncertainty connected with the choice of x* and xl 
is minimized. We note that for ak{x*) {k — 1, 2, 3, • • •) the integration is necessary 
and 9{x) and 5 functions should be kept in the integrand in general. 

As an example let us calculate the small m expansion of f{m). The coefficient ak 
is given at order L by. 



By substituting x* at L = 17 into ak, we have the following satisfactory approximant 
of /(m) = 1 — m + — • • •; 

0.995251 - 0.970003m + 0.9027W - 0.78284m^ + 0.62233m^ + • • • . (19) 



3 Application to the effective potential 

Having prepared basic analysis, we turn to a model field theory which is of our 
main interest. Consider the Gross-Neveu model at the leading order of large N 
expansion^. The Lagrangian is given within dimensional regularization^ sX D — A—2e 

by 

\ g - 

^ V -/v 

£ct = Aa-B^a', (20) 

where 

A=-^^EmLl^., B=fi i = i-^ + log(4.). (21) 

Zn e Zn e e e 



Here MS scheme^ was used for the subtraction. It is well known that the model gen- 
erates the dynamical fermion mass, nidyn — A, where A denotes the renormalization 



group invariant scale in MS scheme. 

At the leading order oil/N expansion, the effective potential is given by the sum 
of diagrams shown in Fig. 2. The straightforward calculation gives 

(22) 

We note that although the naive power counting with respect to leads that the 
contribution with many cr-legs corresponds to higher order in 1/A^, they must be 
included since the vacuum value of a is of order -\/]V- 

The scries (22) converges only when \ga jm\ < 1 and hence the small m behavior 
relevant to the dynamical mass generation is not known from (22). However, Heavi- 



side transformation enlarges the convergence radius and enables us to study the large 
X behavior of the corresponding Heaviside function, V{a,x), as we can see below. 

To obtain x) we need to know the transform of m'^( A; = 0, 1,2, • • •), mlogm, logm 
and 1/m''. Here the following formula is basic. 



n dQ{l/x) 

mil[m) j-^ — -. (23) 

ox 



For example from (1) we have 



log(m) ^ (-7 - log{x))e{x). (24) 



The use of (23) on (24) then leads to 



n 1 

mlogm — > 9{x) — {'J + logx)S{x), (25) 

X 

1 2 / 

m^logm — > —^0(x) 5(x) — H + logx)S (x). (26) 

x^ X 

k 



The transformation of m is easily obtained from 

1 ^ e{x), (27) 

as 

^kH^{k-i)^^^ = 1^2,3, ■■■). (28) 

The 5 functions are needed when one carries out Laplace integrals for the Heaviside 
functions. This is because the 5 function terms cancel out the divergences coming 
from the first terms of (25) and (26), for example. Since the integration over x is 
however not necessary as long as the approximation in the massless limit is concerned, 
we omit, for a while, 5 functions and set 6{x) = 1 in the transformed functions. 
Now, using the results, (24), (25), (26), (27), (28) and Hm-^ = x^/k\, we find 

Note that, due to the creation of l/Zc! in H[l/m''], the series converges for any large 
X. Therefore the large x behavior of V can be easily accessed by increasing the order 
of expansion. This is one of the advantages of V over V. 



We turn to the approximation of the massless effective potential by perturbative 
series at order L, Vl. At L-th order, we have just first L + 1 terms of (29) and find 

X) = ^ - ^ + log A. - , + 1) - t (30) 

The input x* will be determined as in the previous section. Actually the break down 
appears as the domination of the last term in Vl which shows up as its unlimited 
behavior for large x. This can be seen in Fig. 3. And before the domination the 
function experiences a stationary behavior for odd L and large even L. We find that 
the plateau region represents the end of the reliable perturbative regime and thus we 
fix X* by the equation, 

^%^ = 0. (31) 

dx* ^ ' 

If there are several solutions we should input the largest one into x* due to the 
obvious reason. 

Now, the condition (31) gives x* as constant/ ga for odd L . For odd L, the 
substitution of the solution into Vl gives the optimized potential Kpt(c)- For example 
for L = 3, we have the solution, gax* — 1.59607, and this gives the optimized 
potential, 

2 2 

Kpt = ^ (log ^ - 0.373264) . (32) 

The dynamical mass is given from Vopt in the standard way. Note that, since x must 
be positive (see section 2), the region of V^pt thus approximated is restricted to the 
positive (J. In the following, we summarize the result of the approximate calculation 
of the dynamical mass for L = 3, 5, 7, 9, 11; 

L = 3, mdyn/^ = 0.880966 (at gax* = 1.59607), 
L = 5, mdyn/A = 0.97760 (at ^ax* = 2.18061), 

L = 7, mdyn/A = 0.99401 (at gax* = 2.75900), 

tPor even L, there is no solution for (31). In these cases, however, the use of d'^VL{cr, x)/dx'^ = 
gives a solution by which we still have good approximation of rridyn at i = 4, 6, 8, • • •. 
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L = 9, m^ynl^ = 0.99809 (at ^(tx* = 3.33355), 
L = ll, rridyn/^^ = 0.999326 (at ^crx* = 3.90545). 



(33) 



The above result is quite good. Thus, via Heaviside transform approach, the dynam- 
ical mass is approximated only from perturbative information. 

The small mass expansion can be also obtained. Our task is just to substitute the 
solution of (31) into the approximate coefficients. 



ak{x*) 



dx{—x) 



dx 



(34) 



To perform integration, we need the full form of Vl including the 9 and 5 functions. 
The full form is given by 

1 



Vl{x) 



1 

9^ 
27r 



^eix) - -5{x) - (7 + logxpi + 1/2)6' (x) 

h+iogxA-i)«w-|;t|^<)W. 



71 



-6{x) + (7 + logxijL)6{x) 
(35) 



From (35) is given at L-th order as 



TT 



1 

a2 — — 

TT 



Oik 



-- + 1 - 7 + ^og{ga/fi) -logX+Y: 

L / 1 \n Vl' 

\og{ga|^^)-\ogX + 2 + Y. 

n=l 

EriTZ^^V^ (^>2), 



n\n 



where X — x*ga. At L = 11, for example, we have 



(36) 



V{a, m) 



2„2 



27r 

3 



(log^ - 0.499326)+m— (log c/cr- 0.00139409) + — -(log c/a + 1.00595) 



.9o_ 

TT 



W? 1 



2 TT 



+ 



-0.973714 



m 



-0.878954 + 0(m^). 



Qga ' 24c/2(72 

This is quite accurate because the exact result reads 



(37) 



V{a,m) = ^i— (log^--)+-m5(7log^+— (log^+1) 



6^(7 24^V2 



+0{rrt'). 
(38) 



Before closing this section, we prove that our approximants for the massless po- 
tential converges to the exact result in the L ^ 00 limit. That is, limL_»oo Vl((7, x*) — 
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V{a, 0). From (30) we find that lirriL^oo OVl/Qx can be easily summed up and, using 
lim^^oo Vl = V (given by (29)), 

dV{a, x) 1 



exp 



—gax 



(39) 



dx nx^ 

The perturbative truncation of (39) is given by expanding exp[—gax] to relevant 
orders. Since p = cxo for the series expansion of (39), the solution of the truncated 
version of (39) = approaches to oo in the L — > oo limit. More precisely we find the 
scahng of the solution for large L, 

gax* ~ ^L. (40) 
Now consider the reminder Rl, defined by 

n=L+l V / 

Since Vl + Rl = Kx. and p = oo for V^, it is sufficient to show that 

RL{a,x*)^0 (L^oo). (42) 

Note here that x* depends on L and behaves at large L as (40). Now, using the 

Stirling's formula and (40), we have 

,A , 1 (egax*/L)^+'^ 9(gay , 7/0, ,^,r ^ 

\Rl\< 1 \r < -^=^ L-^/^(e/3)^ ^ (L^oo), 

' ' x*V^2^ l-egax*/L V2^(3/e - 1) ^' ' ^ ' 

(43) 

which proves the convergence. 



4 Discussion 

One reason of the success of our approximate calculation is that the transformed 
function has infinite radius of convergence. The other reason is that, as the order 
increases, Vl(c, a;) quickly approaches to the value at x = 00 for fixed a. If one uses 
the closed form, 
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one finds the reason by expanding (44) for large gax, 



V{a,x) 




(log-r — -)+ — ^exp -gax 



gax 



f(-3)(^)V(-3)(-4)(J-f+0((^)^)y 
gax gax gax J 



(45) 



By contrast, the original function has the power- like expansion as shown in (38). 
Thus it is obvious that the approximation of the massless potential is more convenient 
in Vl since it approaches to the "massless" value much faster than Vl- The reason 
behind why the transformed function behaves so good is not known to us. 

We have shown that the perturbative series at finite orders produces the approx- 
imate massless effective potential and the dynamical mass. The deformation of the 
effective potential was made by Heaviside transform with respect to the mass and 
the stationarity prescription to fix the input value x* has found to work good. It was 
also shown that our scheme is capable of approximating the small m behavior. Thus 
the Heaviside transform drastically improves the status of the perturbative approxi- 
mation of physical quantities. We are under the study of full approximation method 
by which the general case where the explicit mass is not small can be treated. The 
result of investigation will be reported elsewhere. 
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Figure Captions 
Figure 1 Model series fiix) is shown for L = 1, 6, 11. 

Figure 2 The Feynman diagrams contributing to the effective potential at the lead- 
ing order oil/N expansion. 

Figure 3 The Heaviside functions of the effective potential with fixed cr at L = 3, 6 
and 9. For the sake of the simphcity, we have set that ga — 1 and A = 1. 
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